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I. INTRODUCTION 



Gravitational wave detection is a very topical subject and great efforts are underway 
in that direction. There are five cryogenic resonant-bar detectors of gravitational waves 
in operation, sensitive ground-based laser- interferometers will be operational soon (LIGO, 
GEO600 and VIRGO) 

m 

meanwhile TAMA300 is collecting data at their level of 
sensitivity and the space-based Laser Interferometer Space Antenna (LISA) is planned to 
be launched around the year 2011 4]. Plans for advanced ground-based interferometers are 
also being developed. 

The most serious hurdle in the detection arises from the very weak interaction of grav- 
itational waves with matter and other fields. Only at energies of order of the Planck scale 
(~ lO^^GeV) will gravitational interaction become as strong as the electromagnetic one. 
Paradoxically this fact makes their detection so attractive, as it would give us informa- 
tion about the epoch at which gravitational waves were decoupled from the cosmological 
dynamics, i.e., the Planck time. 

Gravitational waves can be produced by local sources, such as coalescing stellar-mass 
black holes, compact binary stars and supernovae explosions, or can have a cosmic origin, 
as the decay of cosmic strings and the amplification of zero-point fluctuations due to the 
expansion of the Universe. Waves of the latter type are usually called relic gravitational 
waves (RGWs) and should form an isotropic stochastic background somehow similar to the 
CMBR, but without a thermal distribution 5]. As is well known, RGWs are a unavoidable 
consequence of general relativity and quantum field theory in curved space-time. 

From the point of view of cosmologists obtaining data of RGWs spectrum will be ex- 
tremely interesting as it would make possible to reconstruct the scale factor of the Universe. 
Different cosmic stages of expansion with different equations of estate for the matter con- 
tent will produce a well determined power spectrum of RGWs. Even without data of the 
spectrum of RGWs, we can still get useful information by analyzing their influence in well- 
known processes as primordial nucleosynthesis and the anisotropy they would induce in the 
CMBR. It is possible to test the validity of cosmological models (or to restrict them) by 
comparing the theoretical spectrum of RGWs produced with the maximum values allowed 
by the CMBR anisotropy data and primordial nucleosynthesis. 

The RGWs production has been usually linked to inflation, though inflation is not the 
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Ill ] and a general expression 



/laia jll|. 



Other cosmological 



only process leading to a RGWs spectrum j5|,|6|. By contrast the recently proposed scenario 
of Khoury et al. [7] does not lead to a spectrum of RGWs. The spectrum in an expanding 
universe with a de Sitter inflationary era followed by a radiation dominated era and finally by 
a dust dominated era was calculated by several authors 8|, 
for the creation of RGWs in a multistage model was derived by Maia 
models can be used to obtain in each case a spectrum of RGWs, e.g. an initial inflationary 
stage different from the de Sitter one a model with an intermediate additional era of 
stiff matter between a power law inflation and the radiation era [l^ . 

The aim of this paper is to calculate the power spectrum of RGWs in a spatially flat 
Friedmann-Robertson- Walker universe that begins with de Sitter era, followed by an era 
dominated by a mixture of mini black holes (MBHs) and radiation, then a radiation domi- 
nated era (after the MBHs evaporated) and finally a dust dominated era. We will also see 
what constraints on the free parameters of the cosmological model can be drawn from the 
CMBR. 

As it turns out, the era dominated by a mixture of MBHs and radiation leads to a 
significantly reduced power spectrum of RGWs at present time. This may have implications 
in case LISA does not detect a spectrum. 

Except occasionally, we will use units for which h = c = kB = 1- 

II. RGW SPECTRUM IN AN EXPANDING UNIVERSE 

A. Basics of RGWs amplification 

We consider a fiat FRW universe with line element 

ds^ = -dt" + aitf [dr^ + r^d^^] = a{r]f[-dr]^ + dr"" + r'c/fi'], 

where t and rj are, respectively, the cosmic and conformal time {a{ri)dri = dt). 
By slightly perturbing the metric {g^j = Qij + hij, \hij\ ^ \9ij\j hJ — 

0,1,2,3) the 

perturbed Einstein equations follow. To first order, the transverse-traceless tensor solution 
which represents sourceless weak GWs can be expressed as jl^ |l^ 
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/^S(^,x) = ^G„,(k,x), (1) 



where Greek indices run from 1 to 3, and k is the commoving wave vector. The functions 
GQ,^(k, x) and fi{ri) satisfy the equations 

G^:^ = -k'G^a, Gl^p = = 0, (2) 

/(^) + ik' - ^) KV) = 0, (3) 



a{r]) 

where the prime indicates derivative with respect conformal time and k = |k| is the constant 
wave number related to the physical wavelength and frequency hj k = 27ra/A = 27ra/ = a u. 

The functions are combinations of exp(±ikx) which contain the two possible polar- 
izations of the wave, compatibles with the conditions 

The equation Q can be interpreted as an oscillator parametrically excited by the term 
a" /a. When fc^ ^ i.e., for high frequency waves, expression (jS)) becomes the equation of 
a harmonic oscillator whose solution is a free wave. The amplitude of h^api^i x) decrease 
adiabatically as in an expanding universe. In the opposite regime, when k"^ -C the 
solution to (jni) is a lineal combination of oc aij]) and /i2 oc aij]) J drj a~^. In an expanding 
universe /xi grows faster than fi2 and will soon dominate. Accordingly, the amplitude of 
h^^p{ri,yi) will remain constant so long as the condition fc^ -C ^ is satisfied. When it is 
no longer satisfied, the wave will have an amplitude greater than it would in the adiabatic 
behavior. This phenomenon is known as "superadiabatic" or "parametric amplification" of 
gravitational waves [sl lisj]. 

For power law expansion a oc 77' (/ = —1, 1,2 for inflationary, radiation dominated and 
dust dominated universes, respectively) the general solution to equation is 

/i(r/) = (A;r/)5 J,_i (A;r/) + ir2J_(i_i)(A;r/)) , 

where Ji_^i{kri) , J_^i_i-^{kri) are Bessel functions of the first kind and K12 are integration 
constants. 



B. Spectrum 



Here we succinctly recall the work of Grishchuk 



to obtain the ROW spectrum in a 



universe that begins with a de Sitter stage, followed by a radiation dominated era and a dust 
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era till the present time. Transitions between successive eras are assumed instantaneous. 
The scale factor is 

f 

-ik^ (-00 < < r^i < 0), 

H^iv-'^Vi) {vi<v<m), (4) 



a{ri) 



1 {v+V2-iVl)^ 



{r]2<V < Vo) 



where the subindexes 1, 2 correspond to the sudden transitions from inflation to radiation 
era and from radiation to dust era, Hi represents the Hubble factor at the end of the 
inflationary era and the subindex zero indicates the present time. 
The exact solution to the equation (jH)) for each era specifies to 



cos{kri + — — sin^kr] + 0/) 



Ci 

CRsm{kr]R 



C 



D 



COs{kT]D + 4>d) 



kr]D 



sin{kr]D + 4>d) 



(inflationary era) 
(radiation era) 
(dust era), 



(5) 
(6) 
(7) 



where Ci^r^d, 4>i,R,d are constants of integration, rjR = r] — 2rji and rjr) = rj + rj2 — 4r7i . 

It is possible to express Cr, (j)R and Co, (pD in terms of Cj, 0/ and Cr, (j)R respectively as 
fi{ri) must be continuous at the transition times rj = rji and rj = ri2. Averaging the solution 
over the initial phase 0/ the amplification factor is found to be 



R{k) 



C 



M 



Ci 



,-2 



1 {k » -l/r^i), 
(-l/r^i > A; > l/(r/z)2)), 
k-^ (1/(^D2) > k). 



It is usual to consider that the initial vacuum spectrum of RGWs hi{k) is proportional to 
k j^. The rationale behind that is the following. One assimilates the quantum zero-point 
fluctuations of vacuum with classical waves of certain amplitudes and arbitrary phases; 
consequently it is permissible to equalize hw/2 with the energy density of the gravitational 
waves, c^h^ / {GX^) times A^. Thereby the initial vacuum spectrum of RGWs is given by 
h{k) oc A; B, Q. 

Accordingly, the spectrum of RGWs in the dust era will be 



hf{k) = R{k)hi{k) 



k (A;>-l/r7i), 
(-l/r^i > A; > l/(r/D2)), 

{1/{VD2) » k). 



(8) 
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C. Bogolubov coefficients 



The classical amplification mechanism may be seen as a spontaneous par ticle creation 
by the gravitational field acting on the initial quantum vacuum jiol 11, Q, ^|. This is 
possible as the GW equation is not conformally invariant when a" /a 7^ even in an isotropic 
background ISl. 

□ 

The GW equation may be interpreted as the massless Klein-Gordon equation [10|. Its 
solution can be written as 

/i,,(r/,x) = j (k,x) + (k,a;)) d% (9) 



h^^ (k,a;j- ^e,,(kj ^^^^ e 

where v4(fc), y4|^^ are the annihilation and creation operators, respectively, ejj(k) contains 
the two possible polarizations of the wave and fi{k){v) is a solution to the equation (jH)), as 
in the classical approach, but with the additional condition 

/^{fc)Ai{fc) - f^lk)f^[k) = h (10) 
which comes from the commutation relations of the operators v4(fc), y4|^^ and the definition 



of the field /iij(r/,x) [iall2| 



The solution Q can be expressed in terms of another family of orthogonal modes fi(^k){v) 
as, in contrast with the Mikowskian space-time, in a curved space-time there is no privileged 
family. The two families of modes are related by a Bogolubov transformation 

where a, (3 are the Bogolubov coefficients. The number of particles of the family A^^,, 
contained in the vacuum state of the modes Ji{k) is given by 

(0|iV,|0) = (0|A{,)A(,)|0) = |/3r. 

It is necessary to ascertain which of the solutions to correspond to real particles. An 
approach to that problem is known as the adiabatic vacuum approximation 16]. Basically, 
it assumes that the curved space-time is asymptotically Mikowskian in the limit k ^ 00. 
In that limit, the creation-destruction operators of each family exactly correspond to those 
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associated to real particles and consequently a = 1 and f3 = as the two families of modes 
represent the same vacuum state in the Mikowskian space-time. This argument is identical 
to that made in section II. A, solutions to equation (jS)) with k ^ a" /a are free waves and 
they do not experienj;^ amplification. 

Following Allen |10| and Maia we will now evaluate the number of RGWs created 
from the initial vacuum state in an expanding universe. In the three-stage cosmological 
model we are considering with the scale factor given by (0)), the initial state is the vacuum 
associated with the modes of the inflationary stage fii{ri), which are a solution to the equation 
dni) compatible with the condition (fTUjl . Taking into account the shape of the scale factor at 
this era the modes are 

= {^/2)e'^^k-'/^x'/^H%^{x), (11) 



where x = kr] and ipi is an arbitrary constant phase and Rfl.Jx) is the Hankel function of 



-3/2^ 

-3/2'' 



order —3/2. The proper modes of the radiation era are 

/x« = iV^/2)e^^-k-'/'xfH[%{xn), (12) 

where xr = k{ri — 2?7i) and ipn is again a constant phase. 
The two families of modes are related by 

/i/(r7) = aifiR{r]) + f3ifi*ji{ri) . (13) 

From the continuity of fi{ri) at the transition time rji we obtain 

"^ = -^ + ^^ + ^(4?' = ^^^^ 

where we have neglected an irrelevant phase. The modes with frequency at the transition 
time / = 2TTk/a{rii) larger than the characteristic time scale of transition are exponentially 
suppressed. The characteristic time scale is usually taken to be the inverse of the Hubble 
factor at the transition, H^^ in this case. The coefficients will be ai = 1 and Pi = for 
RGWs with k > 27raiifi and jH)) when k < 2'KaiH. 
In the dust era {rj > ri2) the solution for the modes is 

fiD = {V^/2)e'^^k-'/'x]i'Hj^%{xn), (15) 
where xd = k{ri + ri2 — 4?7i) and is related to the radiation ones by 



Similarly one obtains 

^2 = I 1 



1 



/32 



8{k{7]2-2r],)) 



2 ' 



2k{r]2-2r],) 8 {k{r]2 - 2r],))\ 
for k < 2iTa{ri2)H2 and 0^2 = 1, = for k > 2-na{rj2)H2 where H2 is the Hubble factor 
evaluated at the transition 772- 

In order to relate the modes of the inflationary era to the modes of the dust era, we make 



IG 



For k > 2TTaiHi, we find that 



use of the total Bogolubov coefficients aTr2 and PTr2 
otTr2 = 1, PTr2 = 0; in the range 27raiifi > k > 27ra2-ff2, the coefficients are aTr2 = ai and 
PTr2 = Pi, and finally for k < 2Txa{ri2)H2 we obtain that 

2(r;2 -2r7i) +r/i 



/5' 



'Tr2 



8A;V(r^2-2r/i)2- 

Thus the number of RGWs at the present time 770 created from the initial vacuum state 
is (A^,^) = |/?Tr2|^ ~ ^^'^{Vo) for Lu{r]Q) < 2n{a2/ao)H2, uj-^{r]o) for 27i{ai/ao)Hi > uj{r]o) > 
2vr (02/00) -f^2 and zero for uj{riQ) > 27r(ai/ao)i/i, where we have used the present value of the 
frequency, uj^rjo) = k/aQ. 

Assuming that each RGW has an energy 2huj{ri), it is possible to express the energy 
density of RGWs with frequencies in the range [u{ri),u{ri) + du^i])] as 



dpg{ri) = 2huj{r]) 



27r2c3 



diu(r]) 



(17) 



where P{u{ri)) = {oj'^ijf) /•n'^) {N^) denotes the power spectrum. As the energy density is 
a locally defined quantity, pg loses its meaning for metric perturbations with wave length 
A = 2TT/uj{r]) larger than the Hubble radius H~^{ri). The present power spectrum of RGWs 
predicted by this model is 



P(uj) ^ ( 



(18) 



{uj{r]o) > 27r{ai/ao)Hi) , 
uj~^{r]o) (27r(a2/ao)i^2 < ^^(^70) < 27r(ai/ao)i^i) , 
uj~'^{r]o) {2-nHQ < uj{t]o) < 277(02/00)^2) • 

We can compare these results with the classical ones of the previous subsection. The energy 
density of created RGWs in the classical approach is defined from the spectrum as Q 



Pg 



hi{k)kdk. 



(19) 



27raQ H(i 



From equations (jl7j) and (jl9|) it follows 

P{k) ~ kh){k:). 
8 



Thus both descriptions agree with regard to the dependence on hj and the ranges of the 
wave number of the spectrum of RGWs, as for a power law scale factor the Hubble factor is 
Hiv) = l{ria{ri))-\ 

III. RGWS IN A FRW UNIVERSE WITH AN ERA OF MINI BLACK HOLES 
AND RADIATION 

In this section we evaluate the spectrum of RGWs using the method of the Bogolubov 
coefficients in a more general model than the conventional three-stage model of the previous 
section. 

As is well known, mini black holes can be created by quantum tunneling from the hot 
radiation some cosmological consequences of this effect have been studied 
shall assume these mini black holes (MBHs) are created right after the inflationary period 
(once the reheating is accomplished) and coexist with the radiation until they evaporate. 
During that era the total energy density can be approximated hj p = pbh + Pr and the 
total pressure is 

P = Pbh +Pr = {i- (20) 

where 1 < 7 < 4/3. MBHs follows the dust equation of state pbh = as they can be 
considered non-relativistic matter. If the density of MBHs is large enough to dominate 
the expansion of the Universe, then 7 ~ 1. In the opposite case, the Universe expansion 
is dominated by the radiation, 7 ~ 4/3 and the Universe undergoes the three stages of 
the previous section. From the Einstein equations and (j^UI) one finds 0(77) oc r]\ where 
/ = 27(37- 2) (^ 1 < / < 2) during the "MBHs+rad" era. The MBHs eventually evaporate 
in relativistic particles after some time span whose duration is model dependent. 

A. Power Spectrum 

Assuming that the expansion of the Universe is ab initio dominated by the vacuum energy 
of some scalar field (the infiaton), then dominated by the mixture of MBHs and radiation, 
later radiation dominated (after the evaporation of the MBHs), and finally dust dominated 
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up to the present time, the scale factor of each era is 



a{ri) 



'Wrn {—oo < T] < rji < 0) , (inflation) 



il+T 



m 



Vr 



{rii<ri< r]2), 
{ri2<r] < r]s), 



;"MBHs+rad" era) 
(radiation era) 



^ 2Hi(-r?i)'+^-'?D3 [^^]^ (^3 < ^ < Vo), (dust era) 



where tjbh = V - + l)'?!^ Vr = V + 
r]R2 = [^72 - (1 + l)^i] /I and r]D3 = 2 



(1-0 n+i) , , o(i-0 oC+i) 

. As in the previous section, the 



I (1-0 



transitions betwen stages are assumed to be instantaneous. 

The shape of ^{rj) can be found by solving the equation Q in each era. For the de Sitter 
era ^{rf) is given by (fTTjl as above. For the "MBHs+rad " era the solution of (jH)) is 

where xbh = krisn, and it is related to the modes of inflation by the Bogolubov trans- 
formation (|T3|l with fiR replaced by ^bh- By evaluating the Bogolubov coefficients at the 
transition we obtain 



2o« 



r'2 



when k < 2'KaiHi and a\ = 1, /5{ = when k > 27iaiHi. 

The solution for the radiation era is again |T2|) with = krjn- The coefficients that 
relate (jl2j) with ^bh are 



On 



/32 



1 lTxlxa2 



XR2 



Hf\{lxR2)-H\'MlxB2 



(2) 



1 jTdXR2 



^ +l]H^I'\{lXR2)-Hf^,{lXR2] 



XR2 



when k < 27ra2-ff2 and a2 = 1, P2 — ^ when k > 27ra2-ff2- 

The modes of the dust era are given by ()15|) with xd = krjo and are related to the modes 
of the radiation era by the coefficients 



"3 



1 + 



/?3 



^•^03 



when k < 2TTa^H3 and = 1, jS^ = when k > 2Tra^H^ 
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In order to evaluate the spectrum of RGWs the total Bogolubov coefficients are needed. 
The coefficients relating the initial vacuum state with the modes of the radiation state can 
be evaluated with the help of the relationships 

and one obtains 



1,0 (A; > 2-KaiHi) , 
a\, [3[ {2'KaiHi > k> 2^2/^2) , (21) 



Finally, the total coefficients that relate the inflationary modes with the modes of the dust 
era evaluated from 

are found to be 

«T.3, f^TrS ^ <^ Z-'+^(2^^-3Z+l) ^-(2^+1) ^ rj . ^^^^ 

We are now in position to calculate the current spectrum of RGWs. Taking into account 
that 

-r/i = {aiHi)-\ riR2 = {ai/ aif/\aiHi)-^ , r]D3 = 2{a^/ a2){a2/ aif'\aiHi)~^ , 
and uj = k/ao, the power spectrum P{uj) can be written as 
P ~ (^uj> 2n^Hi 



P ~ ( ^) i72«+2 (2Z-1) (27r^H,>u> 2n^H2] , 

p ^ l^-i2Pjl^l? UV' UV'"' ^4.^-(4.-3) U^a^H2 > > 2n^H- 



(23) 



167r2 



Comparing the power of u in (fTHj) and for < 27i^Hi, we conclude that the 
four- stage scenario leads to a higher number of RGWs created at low frequencies than the 
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three-stage scenario. This fact can be explained intuitively with the classical amplification 
approach. In the three-stage scenario the RGWs are parametrically amplified as long as 
k"^ < a" {rf} / a{ri) . For rj = r]i, a" {rf) / a{ri) vanishes and there is no more amplification. On 
the other hand, in the four-stage scenario the RGWs with uj < 27r^/Ji are amplified until 
rji by the same term a"{ri)/a{ri) = 2/rf than in the three stage model and from rji to 772 by 
the term /(/ — l)/rf^fj. Consequently they have a larger amplitude in the radiation era. 

Figure n] shows the spectrum ()23p for / = 2 and / = 1.1. As is apparent the four-stage 
scenario gives rise to a much lower power spectrum than the three-stage scenario assuming 
that in each case the spectrum has the maximum value allowed by the CMBR bound. The 
higher the MBHs contribution to the energy density, the lower the final power spectrum. 

In this four-stage cosmological scenario, the Hubble factor Hij]) decreases monotonically, 
while the energy density of the RGWs for rj > rj^ can be approximated by pg{ri) ~ if^'*'+^(?7) 
thereby it increases with expansion ll|. Obviously this scenario will break down before 
Pg{f]) becomes comparable to the energy density of matter and/or radiation since from that 
moment on the linear approximation on which our approach is based ceases to be valid. 
The parameters of our model will be constrained such that this cannot happen before the 
current era. 



B. Evaluation of the parameters 

At this stage it is expedient to evaluate the parameters occurring in ((221) • The redshift ^, 
relating the present value of the scale factor with the scale factor at the transition radiation- 
dust, may be taken as 10^ ^]. The Hubble factor Hi is connected to the density at the 



infiationary era by 



c 3mp; 2 



where we have restored momentarily the fundamental constants. In any reasonable model 
the energy density at that time must be larger than the nuclear density (~ lO^^erg/cm^) 
and lower than the Planck density (~ 10^^^ erg/cm^) 22^, therefore 

lO^s"^ <Hi< 10^^s-\ (24) 
Using the expression of the scale factor at the 'MBHs+rad' era in terms of the proper 
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a, ( . / + 1 ^ 



time, one obtains 

where r is the time span of the 'MBHs+rad' era which depends on the evaporation his- 
tory of the MBHs. The simplest model assumes that the MBHs are all created in- 
stantaneously at some given time with the same mass and terminate their evaporation 
simultaneously. If the MBHs evaporated freely, their mass would evolve according to 
{dM/dt)j^^^ = — (7*mp^/(3M^), [2^, where g^: is the number of particles for the black hole 
to evaporate into. However, it is natural to assume that the MBHs are surrounded by 
an atmosphere of particles in quasi thermal equilibrium with them. Therefore we have 
\{dM/dt)^^^\ <^ {dM/dt)j^^^ for a comparatively long period whence the span of the 
"MBHs +rad" era is larger than if the MBHs evaporated freely. Other possible model 
consists in considering that MBHs are created according to some nucleation rate during the 
era, instead of being all nucleated instantaneously. These possibilities give us some freedom 
on evaluating r. 

However, the 'MBHs+rad' era span should be longer than the duration of the transition 
at rji (as the transition is assumed instantaneous) in calculating the spectrum of RGWs. To 
evaluate the adiabatic vacuum cutoff for the frequency we have considered that the transition 
between whatever two successive stages has a duration of the same order as the inverse of 
the Hubble factor. This places the additional constraint 

T > (26) 

Finally, ^ can be evaluated from the evolution of the Hubble factor until the present 
time 

1/2 / N / \ 2 



18^-1 



The current value of the Hubble factor is estimated to be 2.24 x 10 s 



and 



The only free parameters considered here are Z, r and Hi, with the restrictions 1 < / < 2, 
and ()24|). The two first free parameters depend on the assumption made on the MBHs, 
although it is possible to obtain rigorous constraints on Hi and r from the density of the 
RGWs. 
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FIG. 1: Spectrum of RGWs in an expanding universe with a 'MBHs+rad' era for certain values of 
I, T and Hi. The spectrum predicted for the three-stage model of the previous section is plotted 
for comparison, / = 1. It is assumed that each spectrum has the maximum value allowed by the 
CMBR anisotropy data at the frecuency u = 27ri?o = 2.24 x 10~^^s~^. In the bottom-right panel 
the power spectrum with Z = 2 is excluded as it yields a CMBR anisotropy larger than the observed. 



C. Restrictions on the "MBHs+rad" era 



It is obvious that Pg cannot be arbitrarily large, in fact the RGWs are seen as linear per- 
turbations of the metric. The linear approximation holds only for Pg{rj) <C p{r]), p{rj) being 
the total energy density of the Universe. Several observational data place constraints on pg. 
The regularity of the pulses of stable millisecond pulsars sets a constraint at frequencies of 
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order 10~^ Hz [26[. Likewise, there is a certain maximum value for pg compatible with the 
primordial nucleosynthesis scenario. But the most severe constraints come from the high 
isotropy degree of the CMBR. We will focus on the latter constraint. Metric perturbations 
with frequencies between 10^^^ and 10^^*^ Hz at the last scattering surface can produce ther- 



mal fluctuations in the CMBR due to the Sachs- Wolfe effect 
cannot exceed the observed value of 6T/T ~ 5 x 10~^. 



. These thermal fluctuations 



28| 



A detailed analysis of the CMBR bound yields j5, 

Qghlo, < 7 X 10-" i-jj (i/o < / < 30 X H,) (28) 

where Qg = fP{f)/po, po = 'icm?piHl/ {S-nh) and Hq = h^Q^ x 100A;m/(s x Mpc) and 
^100 = 0.7. The CMBR bound for the spectrum evaluated at the frequency u = 271 Hq 
reads js^ 

. ^ 4«+4 /_\2/ rr \3/rr\ 4i-l 

1 > (27r)-^'/2-4i ^2/2-3/ + 1^ 



-4/,2-4«/o;2 Q7 I i\2/'^M (O'0\ ( Hi \ f Hi 



«o/ V«37 V3-72 X 10i9s-V 



and consequently 



/(/, Hi,t) = -mm + I (^28.10 + 2 logio ) - (2/ - 2) log^o (l + ^^i^) (29) 

+ (-4/ + 2) logio(0 + 2 logio(2/2 - 3/ + 1) < 0. 

We next consider different values for / and r. 
(i) When / = 1.1, the relation reads 

/(l.l,i/i,r) = -76.71 + 2.20 logio(^^) - 0.2 log^o (^1 + ^//ir^ < 0, (30) 
see Fig. |21 From it we observe that: 

1. For r < r^=i-i = 1.22 x lO^^^s, the condition is satisfied if Hi < 8.13 x lO^S"^ 
and conflicts with ()26|) . which in the most favorable case is Hi = 8.13 x lO'^^s^^ for 
T = rjr^-^. For / = 1.1, there is no compatibility with the observed CMBR anisotropy 
when r < r^^^-^. Thus, this range of r is ruled out. 

2. For r > rjr^'^, one obtains Hi < H^^^'^{t) from the condition (j30j) . H]r^'^{T) is always 
larger than in the range considered, e.g. H^t^-^It = 10-^°s) = 2.40 x lO'^^s"^. 
Taking into account one obtains that the condition (IHUj) is satisfied for r^^ < 
Hi < H[=^-\r). 
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(ii) When / = 2 (expansion dominated by the MBHs) we have 

/(2, H,, r) = -51.89 + 4 log^o - 2 log^o (^1 + ^^i^) < 0, (31) 

see figure El Inspection of (jM]) and Fig. |21 reveals that: 

1. For r < r^=2 = 6.75 x lO^^S, one obtains Hi < lO^^s ^ which is totally incompatible 
with condition ()26|1 . Hi > 1.48 x lO^'^s"^ for r = rj^^ in the most favorable case. Thus, 
the region r < r^"^ is ruled out as predicts an excess of anisotropy in the CMBR. 

2. For r > r^^^, one obtains Hi < H^t'^It) from the condition (jHTj) . H]r'^{T) is always 
larger than 1.48 x 10^'^s"^ for r in the range considered, e.g. H^r^^T = 10~^°s) = 
1.35 X lO^^s"^. Conditions (f^H|l and (jHlj) are both satisfied in this range for r^^ < 
Hi < Hi=\r). 

We may conclude that the condition ()29|1 leads to different allowed ranges for Hi and r for 
each / considered, although their interpretation is rather similar. For r < rj: the condition of 
minimum duration of the 'MBHs+rad' era ()26|) and the upper bound given by the CMBR 
anisotropy are incompatibles. However, for r > these two conditions are compatible for 
<Hi< H[{t). 

IV. CONCLUDING REMARKS 

We have calculated the power spectrum of RGWs in a universe that begins with an 
inflationary phase, followed by a phase dominated by a mixture of MBHs and radiation, then 
a radiation dominated phase (after the MBHs evaporated), and finally a dust dominated 
phase. The spectrum depends just on three free parameters, namely Hi the Hubble factor at 
the transition infiation-'MBHs+rad' era, r, the cosmological time span of the 'MBHs+rad' 
era, and the power /, being a{ri) oc rf the scale factor of the 'MBHs+rad' era with 1 < / < 2. 

The upper bound on the spectrum of RGWs obtained from the CMBR anisotropy places 
severe constraints on Hi and r. For each value of I considered, there is a minimum value 
of r, r^, compatible with the CMBR anisotropy. There is a range of r, r > r^, for which 
< Hi < HI.{t) satisfies the CMBR upper bound. 
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FIG. 2: The left panel depicts /(l.l, i/i, t)?;s. logio i^i for a) r = 10~''°s, b) r = r^=i-i = 1.23 x 
10~3^s, and c) r = lO'^^s. The right panel depicts /(2,iJi,r) for a) r = lO^^o^^ b) r = r^f^ ^ 
6.75 X lO-^^'s, and c) r = 10-^°s. Conditions (HH) and ^ are satisfied for certain ranges of r 
and Hi in each case. 

The four-stage scenario predicts a much lower power spectrum of RGWs than the con- 
ventional three-stage scenario. We may therefore conclude that if LISA fails to detect a 
spectrum at the level expected by the three stage model, rather than signaling than the re- 
cycle model of Khoury et al. [3] should supersede the standard big-bang inflationary model 
it may indicate a MBHs+radiation era between inflation and radiation dominance truly took 
place. Likewise, once the spectrum is successfully measured we will be able to learn from it 
the proportion of MBHs and radiation in the mixture phase. 

In reality the current epoch is not dust dominated as the recent supernovae type la data 
seems to suggest that the Universe is dominated by a dark energy component and cold dark 
matter implying an epoch of accelerated expansion {2^. This implies a new transition in 
the scale factor shape at some time 174 with tj^ < 1]4^ < tjq. The spectrum of RGWs with 
ko < k < is changed. Thus, the constraints on Hi and r obtained from the CMBR bound 
stay unchanged as the RGWs created at the transition dust-accelerated era were not present 
at the last scattering. 
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